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We present the first calculation of inclusive jet production in deep-inelastic scattering with ap¬ 
proximate next-to-next-to-leading order (aNNLO) contributions, obtained from a unified threshold 
resummation formalism. The leading coefficients are computed analytically. We show that the 
aNNLO contributions reduce the theoretical prediction for jet production in deep-inelastic scatter¬ 
ing, improve the description of the final HERA data in particular at high photon virtuality and 
increase the central fit value of the strong coupling constant. 


INTRODUCTION 

The HERA collider, which operated at DESY from 
1992 to 2007, has produced mauy importaut physics re¬ 
sults, first of all perhaps the most precise determiuatious 
to date of the quark aud gluou deusities iu the protou 
from siugle experiments (HI, ZEUS) [HE] and their com¬ 
bined data sets [3]. These data, taken in deep-inelastic 
electron-proton scattering, are complemented by a wealth 
of data from photoproduction at low virtuality of the 
exchanged photon, in particular on jet production 
giving access also to the distributions of partons in the 
photon [^ and to measurements of the strong coupling 
constant |5j. 

Using the data set of the HERA-H phase of the HERA 
collider from 2003-2007 with an integrated luminosity of 
351 pb“^, the HI collaboration have recently published 
final measurements of inclusive jet, dijet and three-jet 
production in deep-inelastic scattering (DIS) |H jS] and 
used them to determine the strong coupling constant (at 
the mass Mz of the Z-boson) to be 

as{Mz) = 0.1185 ± 0.0016(exp.) ± 0.0040 (th.), (1) 

taking into account absolute double-differential inclusive 
jet, dijet and three-jet cross section data as functions of 
and the jet transverse momentum pp- A more precise 
value was obtained from normalized jet cross sections, 
yielding 

a,{Mz) = 0.1165 ± 0.0008(exp.) ± 0.0038 (th.). (2) 

Unsatisfactorily, only the value obtained by unrenormal¬ 
ized results is in agreement with the current world aver¬ 
age of as{Mz) = 0.1185±0.0006 [S|. The latter uses only 
observables that are known to next-to-next-to-leading or¬ 
der (NNLO) of perturbative QCD, while the analysis of 
the HI collaboration was done in next-to-leading order 
(NLO) accuracy. The lack of knowledge of higher-order 
contributions becomes manifest in a bigger theoretical 


uncertainty due to scale variation in Eq. ([^. The abso¬ 
lute double-differential cross section measurement of in¬ 
clusive jets alone led HI to a value of the strong coupling 
constant of 

as{Mz) = 0.1174 ± 0.0022(exp.) ± 0.0050 (th.). (3) 

In this Letter, we compute the inclusive jet production 
DIS cross section for the first time including NNLO con¬ 
tributions, obtained from a unified threshold resumma¬ 
tion formalism m, and extract an approximate NNLO 
(aNNLO) value for the strong coupling constant. Our 
calculations are based on previous work on inclusive jet 
production in deep-inelastic scattering up to NLO m 
They reduce, as we will see, the theoretical prediction 
and increase the central fit value of the strong coupling 
constant, improving the description of the final HERA 
data in particular at high photon virtuality Q^. 

NNLO CONTRIBUTIONS TO JET 
PRODUCTION IN DIS 

The QCD factorization theorem allows to write the 
differential cross section for inclusive jet production in 
neutral-current DIS with high momentum transfer = 
—q^ as a convolution of the partonic cross section dcr.ya 
with the parton densities in the proton p and the flux 
of photons in electrons as 

do- = X / ^yfi/eiy) 

a 

X J dxp fa/p{xp,pF)d<7 Tr., Tf) , (4) 

where we define y = {p- q)/{p- k) with p and k the mo¬ 
menta of the incoming proton and electron, respectively, 
and q the momentum of the exchanged photon. In deep- 
inelastic scattering the highly off-shell photon has no time 
to decay, so resolved photon contributions can safely be 
neglected. 


2 


From a unified threshold resummation formalism a 
master formula can be obtained that permits to compute 
soft and virtual corrections to arbitrary partonic hard 
scattering cross sections m- At NLO it reads 

daab = dcrf(,-^^ [c 3 -Di(^) + C 2 Do{z) + ciS{l - z)] , 

TT 

( 5 ) 

where for just one color-charged parton in the initial state 
we only need the formula for simple color flow. The func¬ 
tions 


Di{z) 


W(l-z) 
1 - z 


( 6 ) 


with decreasing I are the leading and subleading loga¬ 
rithms at partonic threshold (z —> 1 ) in pair-invariant- 
mass kinematics. The NNLO master formula is given in 
the reference cited above, as are the general formulas for 
the coefficients c^. 

We state here the coefficients for the two partonic 
processes that contribute to jet production in DIS. For 
7 *(; —>■ qg, where 7 * represents the off-shell photon, g a 
gluon and q a quark or an anti-quark, we find 

C 3 = Cf — Nc , (7) 


factors Cp = 4/3 and Nc = 3, on the one-loop /3-function 
/3o = (lliVc —2n/)/3 with n/ the number of active quark 
flavors, the Mandelstam variables t and u, the renormal¬ 
ization and factorization scales gpi and gip and the fixed 
large invariant scale that in DIS is equal to 
whereas in photoproduction it was the Mandelstam vari¬ 
able s. The part Ti of ci does not contain any depen¬ 
dence either on the renormalization or the factorization 
scale. It includes the NLO virtual corrections and is not 
predicted by the threshold resummation formalism. If 
available, it can be read off from a full NLO calculation. 
For the case of transverse photon polarization it could 
be found in reference m- For longitudinal polarization 
we took the formula directly from the source code of the 
corresponding program JetViP mun], that calculates 
inclusive jet production in DIS to NLO accuracy. Some 
two-loop quantities appearing in the NNLO master for¬ 
mula were not given explicitly in | 10 j and could not been 
found in the respective literature. As in our previous 
work on jet photoproduction [5], they were therefore ne¬ 
glected. 


COMPARISONS TO HI DATA 


C2=2C,ln|^)+JVcl„(^ 




and Cl = Cl + Ji with 


^i=-^Cp\n 




-ln(^,. 

4 I M2 / 


For the second process 7*3 —>■ qq we find 
C 3 = 2 {Nc — Cp) , 


( 8 ) 


( 9 ) 


( 10 ) 


C 2 = Nc In 



Nc In 
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and 


c 


M _ 

1 ~ 
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( 12 ) 


These coefficients agree with those found in photopro¬ 
duction for massless jets for the direct part [B]. For mas¬ 
sive jets, additional logarithms depending on the jet ra¬ 
dius R appear (e.g. in C 2 ) [H], which are however irrel¬ 
evant in the case i? = 1 as in the HI analysis consid¬ 
ered here [7]. The calculation is analogous to the case 
of single-jet production in hadron-hadron collisions m- 
The above coefficients further depend on the QCD color 


The NNLO contributions have been implemented in 
the code JetViP for inclusive jet and dijet production in 
DIS, where the convolution over 2 : was already included 
for NLO initial-state corrections on the proton side. At 
NLO, we use of course our complete calculation and not 
only the logarithmically enhanced terms described above. 
As a numerical check, we have repeated the NLO analy¬ 
sis of inclusive single-differential jet production of the HI 
collaboration, performed with NLOJet-|--|- [16] and pre¬ 
sented in reference laiHi, and found excellent agreement, 
confirming previous successful comparisons of different 
NLO programs for jet production in DIS [mill]- 

The measurement took place during the HERA-H run¬ 
ning period with an integrated luminosity of 351 pb“^. 
The beam energies were 27.6 GeV for electrons or 
positrons and 920 GeV for protons, which gives a center- 
of-mass energy of 319 GeV. The leptonic phase space was 
given by 150 GeV^ < < 15 000 GeV^ and 0.2 < y < 

0.7. The jet phase space was restricted to the rapidity 
interval —1.0 < giab < 2.5, where giab is the pseudora¬ 
pidity of a jet in the HERA lab frame. The cross section 
was measured differentially in the jet transverse momen¬ 
tum pt and the virtuality Q^. Jets were reconstructed 
using the fer-clustering algorithm m in the Breit frame, 
where exlusively electroweak processes can be ruled out 
by demanding a minimum of jet transverse momentum 
(here pp > 7 GeV). In inclusive jet production in DIS, 
an almost identical fit result for as{Mz) was obtained 
with the anti-ZcT algorithm. The jet radius was i? = 1. 
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FIG. 1: Single-differential inclusive jet cross sections as a 
function of photon virtuality (top) and jet transverse 
momentum pr (bottom) in NLO (green dashed lines) and 
aNNLO (blue full lines) with the corresponding scale uncer¬ 
tainty bands, obtained by varying pu and pp simultaneously 
by a factor of two up and down, compared to the final HI 
data (red points, color online). 

The perturbative scales were chosen to be 

+Pt)/ 2 and = Q'^ . (13) 

The perturbative hard-scattering functions were convo¬ 
luted with the MSTW2008 set of parton distribution 
functions in the proton with different fixed as{Mz) val¬ 
ues m- This PDF set especially offers the possibility 
to determine a best-fit as{Mz). The number of active 
flavors was nj = 5, since sea contributions of a heavy 
top quark inside the proton can safely be neglected. Fol¬ 
lowing the HI analysis, we choose the PDF member with 
as{Mz) = 0.118 in all plots shown in this Letter, which 
together with the scale choice from Eq. defines our 
central fit. 

In Fig. we compare our NLO (green dashed lines) 
and aNNLO (blue full lines) results to the experimental 
data of the HI collaboration (red points). The uncer¬ 
tainty bands are obtained by varying both scales simulta¬ 
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FIG. 2: Same as Fig. but normalized to the central NLO 
predictions. 


neously about the central scales up and down by a factor 
of two as it was also done by the HI collaboration. We 
have verified that our NLO results and uncertainty bands 
agree very well with those of the HI analysis. Compar¬ 
ing our calculations to the measurements, we see that 
the data points lie in the error bands of both the NLO 
and the aNNLO calculations. For the distribution 
(top), the scale uncertainty is of similar size at NLO and 
aNNLO at low and is considerably reduced in aNNLO 
at higher as expected. We checked that this tendency 
continues and becomes more pronounced for even higher 
Q^, where unluckily we have no data to compare with. 
For the pp distribution (bottom), the scale uncertainty 
is not reduced from NLO to aNNLO, since analytically 
no logarithms of ratios of pp or pp over px appear (see 
above). While in our previous analysis of jet photopro¬ 
duction the jet transverse momentum px ranged from 
17 to 71 GeV [B], its range is restricted in DIS and this 
analysis to lower values of 7 to 50 GeV, while the photon 
virtuality reaches values up to (122.5 GeV)^. 

In Fig. we show the same comparison normalized to 
the central NLO result. We also depict here the central 
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aNNLO result, which is always smaller than the NLO 
central result by approximately 6%. The central dis¬ 
tribution and even more the central pr distribution agree 
better with the HI data at aNNLO than at NLO, as 
they then lie right within the experimental uncertainties. 
While the central NLO results overestimate the measured 
cross sections, indicating that the value of the strong cou¬ 
pling constant as{Mz) used at this order is too large, 
the central aNNLO results underestimate the data and 
require a slightly larger value of Og. The data are thus 
clearly sensitive to the strong coupling constant and can 
be used for an extraction not only at NLO, but also at 
aNNLO. However, we do not expect a significant reduc¬ 
tion of the theoretical error from the scale uncertainty, 
in particular from the pt distribution. 


DETERMINATION OF as 

To determine the strong coupling constant from these 
comparisons, the theoretical calculations have to be per¬ 
formed with a set of parton densities in the proton 
obtained from global fits assuming different values of 
as{Mz)- For our analysis at aNNLO, we employ the 
latest fits of the MSTW group, which have been ob¬ 
tained with NNLO running of the coupling, evolution of 
the parton densities, deep-inelastic scattering and vector- 
boson production matrix elements |19j . 22 different 

MSTW2008 NNLO members were used, which corre¬ 
spond to values of as{Mz) = 0.107 to 0.127. To com¬ 
pare the aNNLO best-fit as to a corresponding NLO 
as based on our full NLO calculation, we carry out the 
same approach with the NLO MSTW2008 parton distri¬ 
bution functions. These are available for the range of 
as(Mz) = 0.110 to 0.130. 

The strong coupling constant as was determined 
by comparing the theoretical predictions at NLO and 
aNNLO to the experimental measurements by HI and 
then finding the minimum value of the reduced x^- 
We present here our results for the single-differential 
PT-distribution, but have verihed that fitting the Q^- 
distribution yields very similar results. At NLO we find 

a^^°{Mz) = 0.115 ± 0.002(exp.) ± 0.005(th.), (14) 

where the central value is slightly lower than the one ob¬ 
tained by HI from the unnormalized double-differential 
inclusive jet cross section (cf. Eq. @) , but where the to¬ 
tal experimental error and the theoretical error obtained 
from a simultaneous variation of the renormalization and 
factorization scales agree very well. At aNNLO, as gets 
shifted upwards, since we demonstrated above that the 
aNNLO contributions reduce the differential cross sec¬ 
tions compared to NLO. We find 

af^^°{Mz) = 0.122 ± 0.002(exp.) ± 0.013(th.), (15) 


where the central value is now slightly above the world 
average, the experimental error is of course unchanged, 
and the theoretical error is slightly larger, reflecting the 
observation made above that the aNNLO calculation is 
not yet sufficiently stabilized by threshold logarithms at 
these values of pr and . The numerical situation would 
only improve at higher values of where unfortunately 
no experimental data are available. 


CONCLUSIONS 

In conclusion, we have presented here a first calcula¬ 
tion of inclusive-jet production in neutral-current deep- 
inelastic scattering up to NNLO of perturbative QCD. 
Leading and subleading logarithmic contributions were 
extracted from a unified threshold resummation formal¬ 
ism for virtual photon-parton scattering processes and 
shown to agree with those appearing in our full NLO 
calculations. The aNNLO contributions implemented in 
our NLO program improve the description of the final 
HI data on inclusive-jet production in the distribu¬ 
tion and even more in the px distribution, when the world 
average value of as is used and for central scale choices. 
The scale uncertainties are reduced only at the highest 
values of where threshold corrections are most impor¬ 
tant. An aNNLO fit of these data with the MSTW2008 
set of parton densities resulted in a new determination 
of the strong coupling constant at the mass of the Z- 
boson that increased the central ht value from below to 
above the current world average, but did not reduce the 
theoretical error. 
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